In the claims: 

1. An image enhancement method, comprising: 
capturing an image; 

constructing a multi-resolution structure comprising one or more resolution layers; 
processing each resolution layer using an iterative algorithm having at least one 
iteration; 

projecting each processed resolution layer to a subsequent resolution layer; 
up-calling each projected resolution layer to the subsequent resolution layer; and 
using the projected resolution layers to estimate an illumination image. 

2. The method of claim 1, further comprising, for each of one or more iterations: 
calculating a gradient of a penalty functional; and 

computing an optimal line-search step size. 

3. The method of claim 2, wherein the penalty functional is given by: 

F [l] = j(v/| 2 + a(l - sf + /3\V{1 - sf)ixdy ; 

Q 

subject to / > s and (vi, n^j - 0 on dQ; wherein Q is a support of the image, 9Q is an image 

boundary, n is a normal to the image boundary, and oc and P are free non-negative real 
numbers 

4. The method of claim 2, wherein the penalty functional is given by: 
F[l]= ^(V^V/| 2 +a(/-s) + j3w 2 (Vs]Vl- Vs\ 2 }lxdy 

where Wi and w 2 are non-linear functions of the gradient. 

5. The method of claim 1, wherein the iterative algorithm is a Projected Normalized 
Steepest Descent algorithm. 

6. The method of claim 1, wherein the iterative algorithm is a Steepest Descent 
algorithm. 

7. The method of claim 1, wherein a set of constraints comprise that the illumination is 
greater than the image intensity, 

8. The method of claim 1, further comprising applying penalty terms, the penalty terms, 
comprising: 

that the illumination is spatially smooth; 

that the reflectance is maximized; 

that the reflectance is piece-wise smooth. 
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9. The method of claim 1, further comprising: 

computing the reflectance image based on the captured image and the estimated 
illumination image; 

computing a gamma correction factor; 

applying the gamma correction factor to the estimated illumination image; and 
multiplying the gamma-corrected illumination image and the reflectance image, 

thereby 

producing a corrected image. 

10. A system, embodied in a computer-readable medium, for enhancing digital images, 
comprising: 

a log module that receives an input digital image S and computes a logarithm s of the 
input digital image; 

an illumination estimator module that produces an estimate /* of an illumination 
component L of the input digital image S, wherein the estimator module employs a construct 
comprising one or more resolution layers, and an iterative algorithm that processes each of 
the one or more resolution layers; and 

a summing node that sums the logarithm s and a negative of the estimate /* to produce 
an estimate r*of a logarithm of a reflectance component R of the input digital image S, 
wherein a processed resolution layer is used to up-scale a subsequent resolution layer. 

11. The system of claim 10, wherein the iterative algorithm, for each of one or more 
iterations: 

calculates a gradient of a penalty functional; and 
computes an optimal line-search step size. 

12. The method of claim 11, wherein the penalty functional is given by: 



subject to 1 > s and ( V/, n) = 0 on 3fi; wherein Q is a support of the image, dQ is an image 



boundary, n is a normal to the image boundary, and a and (3 are free non-negative real 
numbers. 

13. The system of claim 10, wherein the penalty functional is given by: 



F[l]= J(V/| 2 + «(Z 



sf+/3\v(l-sf)cbccfy 



a 




F[l] = ^ (w, (Vs] V/| 2 + a(l -s) + J3w 2 (Vs}Vl - Vs| 2 \ixdy 



where wi and W2 are non-linear functions of the gradient. 
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1 14. The system of claim 10, wherein the iterative algorithm is a Projected Normalized 

2 Steepest Descent algorithm. 

3 15. The system of claim 10, wherein the iterative algorithm is a Steepest Descent 

4 algorithm. 

5 16. The system of claim 10, wherein each of the one or more resolution layers is projected 

6 onto constraints, and wherein the constraints comprise that the illumination is greater than the 

7 image intensity, L>S; 

8 17. The system of claim 10, further comprising penalty terms, the penalty terms 

9 comprising: 

10 that the illumination is spatially smooth 

1 1 that the reflectance is maximized; and 

12 that the reflectance is piece-wise smooth. 

13 18. The system of claim 10, further comprising: 

^3 14 a module that computes reflectance and illumination images based on the input 

Cy 15 digital image and the estimated illumination image; 

^ 16 a gamma correction module that computes a gamma correction factor and applies the 

CS 17 gamma correction factor to the estimated illumination image; and 

iy 

2 18 a node that multiples the gamma-corrected illumination image and the reflectance 

l n 19 image, thereby producing a corrected digital image. 

n ~ 20 19. A method for enhancing an image S, the image S comprising a reflectance R and an 

21 illumination L, the method comprising: 

rg 22 constructing a multi-resolution image structure having one or more resolution layers; 

r=s 23 processing the resolution layers using an iterative algorithm; 

24 projecting the processed resolution layers onto a set of constraints, the set of 

25 constraints 

26 comprising boundary conditions and that L>S; and 

27 using the projected resolution layers to estimate an illumination image. 

28 20. The method of claim 19, wherein the image S is a RGB domain color image, the 

29 method further comprising: 

30 mapping colors R, G, B of the image S into a luminance/chrominance color space; 

3 1 applying a correction factor to a luminance layer; and 

32 mapping the luminance/chrominance colors back to the RGB domain. 
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